Draft version February 5, 2008 

Preprint typeset using LATpjX style emulateapj v. 11/12/01 



IMPROVED FOKKER-PLANCK EQUATION FOR RESONANCE LINE SCATTERING 

George B. Rybicki 

Harvard-Smithsonian Center for Astrophysics 
60 Garden Street, Cambridge, MA 02138 USA, grybicki@cfa.harvard.edu 
Draft version February 5, 2008 

ABSTRACT 

A new Fokker-Planck equation is developed for treating resonance line scattering, especially relevant 
to the treatment of Lya in the early universe. It is a "corrected" form of the equation of Rybicki & 
\^ , Deir Antonio that now obeys detailed balance, so that the approach to thermal equilibrium is properly 

described. The new equation takes into account the energy changes due to scattering off moving particles, 
the recoil term of Basko, and stimulated scattering. One result is a surprising unification of the equation 
for resonance line scattering and the Kompaneets equation. An improved energy exchange formula due to 
resonance line scattering is derived. This formula is compared to previous formulas of Madau, Meikson, 
& Rees (1997) and Chen & Miralda-Escude (2004). 

Subject headings: cosmology — atomic processes — radiation processes: general — radiative transfer 
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1. INTRODUCTION 

The scattering of radiation in resonance lines plays an important role in many parts of astrophysics, including stellar 
atmospheres, diffuse nebula, and active galactic nuclei. In the early universe, resonance scattering in the hydrogen Lya 
' line has been shown to be an crucial process for understanding the state of the primordial gas during the recombination 
, era (Peebles 1968; Zel'dovich, Kurt, & Sunyaev 1969). It can also can play a critical role in determining the populations 
£f) • of the fine structure levels in the 15 level of hydrogen (Wouthuysen 1952; Field 1958, 1959), and thus affect the formation 
\ of the 21cm line, observations of which promise to be an important source of information about the early universe (Hogan 
\D . & Rees 1979; Madau, Meikson, & Rees 1997; Tozzi, Madau, Meiksin, & Rees 2000). It may also play a role in direct 
O ■ heating of the gas (Madau, Meikson, & Rees 1997; Chen & Miralda-Escude 2004). 

The theoretical description of resonance line scattering is based on a redistribution function R, which gives the probability 
O that an initial photon state will be scattered into some final photon state. The level of description of these photon states 
depends on the problem treated. For this paper we confine our attention to problems where the radiation field is sufficiently 
isotropic that only the angle-averaged form of the redistribution function needs to be used, and also that the radiation 
can be considered unpolarized. The redistribution function then depends only on the intitial v' and final v frequencies of 
^ \ the scattered photon, that is, R{y,v'). 

For Lya and other resonance line scattering in low density media, coherent scattering in the atom's rest frame from 
a natural (Lorentz) profile is appropriate. When one accounts for the Doppler effect due to atoms with a Maxwellian 
velocity distribution, a scattered photon will have its frequency decreased or increased depending on the components of 
jJJ ■ the atomic velocity along the initial and final photon directions. The appropriate redistribution function under these 
& ' conditions was first derived by Henyey (1941) and, following Hummer (1962), is usually denoted Rii{v, v'). For many 
cases in astrophysics, Ru captures the dominant physics for resonance line transfer. 

The solution of line transfer problems using a full redistribution function Ru is not trivial, since all frequencies are 
coupled together. However, there is an important class of problems where an approximate formulation can effectively used, 
namely, when the radiation field is sufficiently smooth on the scale of the Doppler width of the line. Then it is possible to 
derive a Fokker-Planck (F-P) type of transfer equation, where the redistribution is taken into account by a second-order 
differential operator over frequency space. An equation of this type was first given by Unno (1952), and improvements 
were made subsequently by Harrington (1973), Basko (1981), and by Rybicki & Dell' Antonio (1994, hereafter RD). 

Another effect on redistribution is the loss of photon energy during scattering due to the recoil of the atom. This is 
completely analogous to the ordinary Compton effect, but is smaller by the ratio of electron to atomic mass. The recoil 
effect was first discussed by Field (1959) for the case of resonance scattering with zero natural line width (i?i in the 
notation of Hummer 1962). Adams (1971) qualitatively considered recoil for Ru in the context of Lya scattering in 
the neutral hydrogen of the galactic disk. Basko (1981) derived the appropriate generalization of the Rjj including recoil 
redistribution function and showed that it led to a simple additional term in the F-P formulation. Rybicki & Dell' Antonio 
(1994) showed how the Basko term could be included in their F-P equation and how it affected its solutions. 

In all the cases treated by Adams, Basko, and RD the recoil term did not seem to be important, at least for its effect on 
the radiation field itself. However, Madau, Meikson, & Rees (1997) argued that the atomic recoil in Lya scattering could 
be an important heating mechanism for the intergalactic medium before reionization, basing their analysis on a simple 
formula that included recoil but no other effect. This was also investigated by Chen & Miralda-Escude (2004), who found 
a much lower heating rate using the RD F-P equations with Basko term. 

A strong motivation for the present paper is to understand better the difference between the results of Madau, Meikson, 
& Rees (1997) and Chen & Miralda-Escude (2004) regarding energy exchange. It seems clear that Chen & Miralda-Escude 
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(2004) must be at least partially right: if the radiation field were thermal at the same temperature as the gas, then no 
net energy exchange should occur, so that the simple formula of Madau, Meikson, & Rees (1997) cannot be of general 
validity. On the other hand, the formulation of Chen & Miralda-Escude (2004) did not take detailed balance fully into 
account, so the matter deserves further investigation. 

Another motivation for this paper is the issue of the 21cm spin temperature and how it depends on Lya transfer. 
According to the Wouthuysen-Field effect (Wouthuyscn 1952; Field 1958, 1959) the radiation field near Lya comes into 
thermal equilibrium with the kinetic temperature of the gas, and this Lya field then imposes that same spin temperature 
on the 21cm line. Any transfer equation, cither using a full redistribution function or through a F-P approximation, should 
correctly incorporate this thermal requirement, since small errors in calculation could have large effects on predictions of 
21cm radiation. 

It should be obvious that an accurate determination of energy exchange rates and 21cm spin temperatures is critically 
dependent on being able to formulate the interaction of of radiation and matter in a thermodynamically correct way. 
Deguchi & Watson (1985) showed that this will be true if the redistribution functions satisfy certain detailed balance rela- 
tions. However, no known form of Ru satisfies these relations exactly, so Deguchi & Watson (1985), in their investigations 
of the recombination era, found a method of "correcting" the usual Ru, such that its values did not change very much, 
but enough to satify detailed balance. Deguchi & Watson (1985) applied their correction method to the redistribution 
functions directly and solved the equations numerically; they did not consider how this might change any associated F-P 
method. 

A review of the basic equations is given in §2. In §3 we introduce a simple correction scheme, different from that of 
Deguchi & Watson (1985), which can be used to produce a redistribution function satisfying detailed balance, and we use 
this correction scheme to derive a new, corrected F-P equation. This F-P equation incorporates the RD diffusion term, 
the Basko recoil term, plus a new term that describes stimulated scattering and another that accounts for the correct 
phase space factors; acting together, these terms ensure detailed balance and the approach to the proper thermal radiation 
field. In §5 we demonstrate the close analogy between the new corrected F-P equation and the Kompanects equation. In 
§6 we use the new F-P equation to derive an energy exchange rate formula. This will be compared to the the formulas of 
Madau, Meikson, & Rees (1997) and of Chen & Miralda-Escude (2004). 

2. BASIC EQUATIONS 

We shall consider the time dependence of the radiation field in a homogeneous and isotropic medium. This idealized 
problem allows one to develop equations applicable to more complex situations under the isotropic scattering approxima- 
tion, much as is done for the analogous Kompaneets equation. When dealing with scattering problems, it is often easier to 
use intensities based on photon numbers rather than energies, and this is done here. Because of the isotropic assumption, 
the radiation is determined by the mean intensity J(v). 

For simplicity of presentation, we shall write all of our transfer equations for the special case of the time dependence of 
a homogeneous, static medium, very much as in the spirit of the usual Kompaneets equation (Kompaneets 1957; Rybicki 
& Lightman 1979). One should interpret the right hand side of such equations as "emission minus absorption" terms, 
which can be incorporated into other transfer equations as needed. In particular, we do not include a term accounting for 
the expansion of the medium, which, if needed, should be included on the left hand side as part of the transfer operator. 
For cases that involve spatial transfer, one can also use the results here in an "isotropic scattering approximation" by 
appropriate modification of the absorption terms, that is, by adding back the isotropic absorption term and subtracting 
the proper anisotropic absorption term. 

For many applications the process of stimulated scattering is not important, and most of the discussions of resonance 
line scattering have explicitly or implicitly neglected it. In that case the time dependence of the radiation is then governed 
by the transfer equation, 

= J [R(u, v')J{v') - R(v', u)J{v)\ dv' = -<fi{v)J{v) + J R(v, v')J{v') dv' . (1) 

Here R(v, v') is the Ru redistribution function and ip(v) is the Voigt line profile function, related to the redistribution 
function by, 

<p(v) = J R{v',v)dv', (2) 

and which is normalized to unity, 

/'OC 

<p(v)dv = l. (3) 
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The quantity \ is given by 



x = ^ Nl B 12 , (4) 

where vq is the line center frequency and B\2 is the upward Einstein i?-coefficicnt for the resonance transition. 

When stimulated scattering is important, it can be taken into account through additional factors in the transfer equation 
(1 + n), where n(v) — c 2 J(v)/2v 2 is the photon occupation number, 

— 1ST = / + R <?> - t 1 + n{v '^ ^ ") J ^ dv '- ( 5 ) 
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It is important to note that the same redistribution functions appear in both equations (1) and (5), that is, they are 
independent of whether stimulated scattering takes place or not. Therefore, it is possible to determine the properties of 
the redistribution functions for cases without stimulated processes, and then use them in equation (5) when necessary. 

Both equations (1) and (5) conserve photons, as can be proved by integration over v and interchanging v and v' in one 
of the two terms on the right. As such, photon conservation is built into the transfer equations themselves, and does not 
depend on a particular symmetry or other property of the redistribution functions. This is as it should be, since scattering 
does not change photon number. 

The difficulty of solving transfer problems using the full redistribution functions has led to some simplified formulations 
that are adequate under special circumstances. One important fact is that the redistribution function is significantly 
nonzero only for frequencies v and v' that do not differ by many Doppler widths. In many cases of interest the radiation 
field in the neighborhood of line center is sufficiently smooth such for a given value of v, the relevant values of J(y') can be 
represented by a Taylor series up to second order in the difference (V — v). This makes it possible to replace the frequency 
integrations on the right hand side of the transfer equation with an equivalent second-order differential operator. The 
resultant transfer equation is known as a Fokker-Planck (F-P) equation. 

There have been a number of Fokker-Planck formulations for the Ru redistribution problem given in the literature. 
They fall into two categories: first there are ones that are based on the standard Ru redistribution function (Henyey 
1941; Hummer 1962), which consider the spreading of radiation on scattering but neglect recoil. These formulations give 
an expression for the integral of the redistribution function and the intensity by means of a differential operator, 

dJ{u)' 
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where R = Ru- Here, 
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is the usual Doppler width and D(v) is the "diffusion" coefficient, essentially a second moment of the redistribution 
function. This implies the F-P transfer equation (without stimulated processes), 
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The different F-P formulations differ in what they take for D(v). The first such formulation by Unno (1955) (see also 
Harrington 1973) took 

which is an asymptotic result for the second moment, strictly applicable only in the wings of the line. Basko used a 
non-singular form of (9), 

DM = MWW ' (10) 

equal to the natural line profile, but did not justify this choice. Rybicki & Dell' Antonio (1994) rigorously caculated the 
second moment of Ru using the full Ru redistribution function. Their fundamental result is given in their equation (A13) 
for the second moment of the redistribution function including angular dependence. Averaging that result over angle 
depends on whether an isotropic or dipole phase function is assumed. The two results are, in terms of physical frequency, 

D(v) = <p(v) + \av^ <p"(v), (RD isotropic), (11) 

D(v) = ip{v) + ^Ai£ <p"(v), (RD dipole). (12) 

In fact, RD adopted a version without the second derivative, 

D{v)=y{v), (RD adopted). (13) 

Their main rationale for this choice was that either the first or second moment needed to be altered in order to preserve 
particle conservation, and (13) was the simplest way of doing this. Others (Madau, Meikson, & Rees 1997; Chen & 
Miralda-Escude 2004; Hirata 2005) have used the RD approximation with this form for D(v). 

For much of formal developments in this paper, D(v) will be left general, and only later will particular choices be made. 
One should point out that all of the expressions (9)-(13) for D(v) are asymptotically equivalent in the line wings, so 
it only near the line core that differences occur. A desire for a non-singular result argues against (9). Although RD is 
somewhat ambiguous as to the correct form, it seems clear that the correct cutoff for the asymptotic wing region should 
occur on the scale of the Doppler width; this argues against Basko's choice (10), which implies this scale is the natural 
width. In §5, the choice for D(v) will be discussed further. 

The second category of F-P equations take recoil into account. This was done first by Basko (1981), who gave the 
following generalization of the F-P equation (8), 
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in which the new term hJ/kT accounts for the recoil. Basko used the value (10) for D(v). RD introduced a slight variation 
of Basko's generalized equation by using (13) for D(v), which was subsequently employed by CM. This equation with 
recoil will not be used directly in our derivation of a corrected F-P equation, although we shall compare our results with 
it later. Unless indicated to the contrary, the term "uncorrected F-P equation" will refer strictly to equation (8). 

3. DETAILED BALANCE AND "CORRECTED" REDISTRIBUTION 

In order to describe correctly the approach to thermal equilibrium, the redistribution function must satisfy a certain 
symmetry condition known as detailed balance, which is a result of the microscopic reversibility of the scattering process, 
along with the inclusion of appropriate phase space factors and the fact that the scattering particles obey a Maxwellian 
velocity distribution at temperature T. 

Deguchi & Watson (1985, hereafter DW) considered the detailed balance relations satisfied by redistribution functions. 
Their result was stated for the redistribution function using energy based intensities, and they also used opposite ordering 
of arguments v and v' . Accounting for these differences, the DW result translates in our notation to 

R(u', vy e - hv / kT = R(v, u'y 2 e- hv 'l kT . (15) 

In most cases of line transfer, one usually deals with frequencies v and v' very close to the line center frequency and 
thus very close to each other. Then equation (15) can be approximated by the symmetry relation, 

R*{v',v)=R*(v,v'). (16) 

We have added the star to denote redistribution functions satisfying the special symmetry (16) rather than (15). In fact, 
the usual redistribution functions used in line transfer calculations are clearly approximate in that they satisfy (16) rather 
than (15) (e.g., Hummer 1962). 

Approximations such as equation (16) are satisfactory for many applications. However, as described in the Introduction, 
there are cases where one needs a formulation that includes enough physics to ensure the detail balance relation (15). 

Let us consider how one may modify a simple redistribution function obeying (16) so that it satisfies (15) instead. 
Deguchi & Watson (1985) did this by choosing the usual expression for v > v' , and obtaining the result for v < v' by the 
detailed balance relation [see discussion in RW after equation (RW22)]. In our notation, this may be written, 

p (,. „>\ _ j R*{vy), v>v', , s 

«dwC, > ~ y R* (v' ,v)(v / v' ) 2 cxv[h(v' - v)/kT], v < v> . (U > 

While this choice clearly will lead to thermal equilibrium in the appropriate limit, it does have the undesirable property 
of creating a discontinuity in slope in the redistribution function at v = u', even when one does not exist in R*. 

In this paper we use a different method of correction. Again starting with an R* that satisfies the approximate detailed 
balance relation (16), we define the "corrected" redistribution function R Cl which is related to the original R* through a 
"diagonal similarity transformation," 

R c {v, v') = U{v)R*{v, v'y- 1 ^'), or R*(u, v') = E7 _1 (i/)i2 c (i/ ) v')U(v'), (18) 

where, 

U{v) = ve-^ hv ' kT . (19) 

It is easily verified that this corrected redistribution function satisfies the detailed balance relation (15). We prefer this 
method to that of DW, because it introduces no new discontinuities of slope at v = v' and is more amenable to the kind 
of analytic manipulations employed in this paper. 

4. CORRECTED FOKKER-PLANCK EQUATION 

The correction proceedure given in equation (18) can be used in any formulation of the redistribution problem. For 
this paper, however, we shall apply it to the Fokker-Planck equation (8). In Appendix A, it is shown that the R(v,v') 
implied by equation (8) satisfies the symmetry condition (16), so that equation (18) may be used to define a corrected 
redistribution function. Using this, we straightforwardly derive our new corrected Fokker-Planck equation that obeys 
detailed balance, namely, 
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dJ ( h 2\ T c 2 J 2 
J 



(20) 



(see Appendix B for the mathematical details of this derivation.) 

We now identify the physical significance of the terms in the inner set of brackets of this corrected F-P equation. The 
first term, dJ/du, is easily understood, since it is the same diffusion term as in the uncorrected equation (8). The last 
term, involving the square of J, clearly describes the stimulated scattering; let us ignore this term for the moment, since 
we want to make comparisons with other F-P equations that have also ignored this process; we shall return to it in §5. 

There are two remaining terms, hJ/kT and —J/2v. The first of these is precisely the Basko recoil term, as in equation 
(14). This may seem remarkable, since the F-P equation in its uncorrected form (8) does not include recoil at all. 
However, there is a simple physical explanation for this seeming coincidence. In the theory of Brownian motion, which 
is also described by a F-P equation, Einstein found that thermal equilibrium required a certain relations between the 
coefficients of "diffusion" and "drift" (other terms for drift are "mobility" or "dynamical friction"). In the present case, 
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the uncorrected equation (8) describes only the "diffusion" part of the process. Imposing detailed balance through our 
correction procedure has automatically generated the proper "drift" term, here given by the recoil term. 

The remaining term —J/2v can be explained as follows. With just the diffusion term and the Basko recoil term, as in 
equation (14), the equilibrium solution (setting the bracket to zero) would give J oc exp(— hv/kT). In fact, Basko (1981), 
(following Field 1959), claimed that his treatment of recoil led to thermal equilibrium because the equilibrium solution 
of his equation was proportional to exp(— hv/kT). However, to be rigorously thermal (without stimulated processes), 
J should be proportional to the Wien law, v 2 exp(— hv/kT), not just exp(— hv/kT). The extra term —2J/v in the F-P 
equation (20) is just what is needed to account for the missing factor of v 2 . In many applications this additional slowly 
varying factor may be unimportant. However, since the goal here is to account for thermal equilibrium as precisely as 
possible, this term should be included as a matter of principle. 



5. RELATION TO THE KOMPANEETS EQUATION 

For further investigation of the corrected equation (20), especially the stimulated scattering process, it is useful to write 
it in terms of the photon occupation number n, 



J_dn 

cx dt 



Av 2 B 



1 d 



2 v 2 dv 
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dv + kf n{1 + n) 



(21) 



This is reminiscent of the Kompaneets equation (Kompaneets 1957; Rybicki & Lightman 1979), which can be written, 
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dn 



kT 1 d 

m e c 2 v 2 dv 



dn 
dv 



kT 



n(l + n) 



(22) 



N c a^c dt 

where N c is the electron density and or is the Thomson cross section. 

Some of the desirable general properties of the new equation follow exactly as in the Kompaneets case. It satifies photon 
conservation exactly owing to the form of the right hand side as 1/v 2 times a frequency derivative. It satisfies detailed 
balance by construction; a result of this is that it has the correct thermal equilibrium solution. We can see this by setting 
the right hand side to zero, which leads to the condition, 



dn h ,„ . 

_ + _ n (l + „) = 0. 

The general solution to this differential equation can be written 

1 

n ■ 



(23) 



(24) 



where fj, is a constant. This is a thermal Bose-Einstein distribution with chemical potential /x, which is the correct 
thermal distribution for photons obeying a conservation law. We note that the general properties of photon conscvation 
and detailed balance arc independent of the particular form of D(v). 

The similarity between the corrected RD equation (21) and the Kompaneets equation (22) can be emphasized further 
by using definition (7) and choosing the RD value (13) for the diffusion coefficient D(v). Then equation (21) can be 
written, 



1 dn 

Vdt 



kT 1 d 

nip-c 2 v 2 dv 



.2.. 2, 



dn 
dv 



h , 



(25) 



where m = m a is the atomic mass and a is the line absorption coefficient, 

a L (v)=xv(v). (26) 

We now make a minor modification of equation (25) by replacing v 2 by v 2 . This is a very slight change, since a(v) is 
significant only in a narrow region around line center, so the distinction between v and vq is slight. Such a replacement 
is not unusual for line problems; it is quite typical when deriving line redistribution functions, and is even found in the 
definition of the Doppler width, equation (7). 

One may also regard the above modification as equivalent to a slightly different diffusion coefficient D(v), and, as 
noted above, the important properties of photon conservation and detailed balance are independent of D(v). With this 
modification, equation (25) becomes 



1 dn 
c~dt 



kT 1 d 

m a c 2 v 2 dv 



4 

v a 



dn 
dv 



+ fcT n(1 + n) 



(27) 



Remarkably, equation (27) is precisely the same as the Kompaneets equation if we take m to be the electron mass m 
and a to be the Thomson scattering coefficient, «t = N c (Tt- Thus, equation (27) provides a unification of the Fokker- 
Planck equation for resonance line scattering and the Kompaneets equation. One only needs to substitute the appropriate 
scattering coefficient and mass of the scatterer. 

It should be noted that this unification is even valid if one chooses to use one of the forms (11) or (12) for D(v) instead 
of (13). In either case, this simply involves replacing a by a + Ca" in equation (27), where C is a constant. This will still 
be completely consistent with the Kompaneets equation, because a for Thomson scattering is independent of frequency. 
However, for the remainder of this paper we shall continue to use the corrected Fokker-Planck equation in the simple form 
(27). 
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6. ENERGY EXCHANGE FORMULA 



Equation (27) can be used to derive an energy exchange formula between the radiation and kinetic energy of the atoms. 
The energy density of the radiation field is 



U 



— / v z ndv. (28) 



Differentiating with respect to time and using equation (27), followed by an integration by parts, yields the rate of change 
of radiation energy density, 



dU _ SnhkT 
dt rriaC 4 



dn h 



du, (29) 



which is the negative of the heating rate for the gas. 

Equation (29) is the general formula for energy exchange within our Fokkcr-Planck formulation. The presence of the 
frequency derivative in the integral makes the formula difficult to interpret, but it may be put into a more revealing form 
by defining at each frequency an appropriate temperature of the radiation field by comparison to a thermal Bose-Einstein 
distribution with temperature T and chemical potential fi, 

" th = e (h,-^/kT _ I for which n > th = -JL n (l + n). (30) 

Given the values at frequency v of an arbitrary radiation field n(y) and its slope n'(u), equation (30) one may define the 
local "radiation temperature" Tr(i/) for an arbitrary radiation field by 

IW „)__*»(L+2!>. (31) 

With this definition, equation (29) becomes 

dU 8irh 2 f , / T 



ot m **j*^« 1+n A l -w)r- (32) 

In order to interpret this equation, let us consider the special case when the radiation field is smooth and slowly varying 
in the neighborhood of the line center at vo, so it and its derivative can be considered constant. Then the delta- function 
character of a^v) implies, 

; n n 1 + nn) 1 - -— 



where the subscript "0" means evaluation at uq. 

This may be put into another form by defining the rate of scattering per unit volume by 

P = 4tt y ol(1 + n)Jdv = ^^n (l + no), (34) 

where we have again used the delta- function character of a^v). Note that stimulated scatterings are included in this 
rate. Then our energy exchange formula takes the simpler form, 



9f - P {^)^{ l -T^)' < 35 > 

(We remark that the general energy exchange formula (32) could be written in this same simplified form, provided one 
uses for the radiation temperature a suitable mean value, frequency-averaged over the line.) 

The factor (1— T/Tro) in the formula (35) shows how the energy exchange between the radiation and the matter depends 
on their relative temperatures. In particular, energy always flows from the hotter to the colder system; no energy exchange 
occurs when the radiation and matter temperatures are the same. These fundamental principles of thermodynamics are 
guaranteed by the detailed balance built into our equations. 

The factor (1 — T/Tro) is particularly important for situations with large optical depth in the resonance line. In 
these cases, the Wouthuysen-Ficld effect (Wouthuysen 1952; Field 1959) becomes relevant. According to this effect, the 
scattering within the line changes the local radiation field in such a way that the radiation temperature within the line 
is brought close to the gas temperature. Therefore the factor (1 — T/Tr ) becomes very small and so greatly reduces 
the flow of energy between the radiation and the gas, to a value much less than one what might estimate from simple 
arguments based on the radiation temperature of any incident radiation. 

We are now in a position to compare our energy exchange formulas to those of Madau, Mcikson, & Rees (1997, hereafter 
MMR) and Chen & Miralda-Escude (2004, hereafter CM). 

6.1. Comparison with MMR 

In deriving their energy exchange formula, MMR assumed (implicitly) the near constancy of the radiation field near 
line center, so it should be compared with our equation (35). Some translation is necessary, since their rates are defined 
per atom instead of per unit volume. Making that adjustment, the MMR formula (their Eq. [30]) is, in our notation, 

■aT = - p \ ZT-JI h »o- (36) 



dt \m 
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The rate P is the same as our rate P, defined in equation (34), except that it omits the stimulated scattering factor 
(1 + n). 

The omission of stimulated scattering is minor, but more importantly the MMR formula does not contain the factor 
(1 — T/Tro), which, as we have seen, describes how the energy exchange depends on the relative temperatures of the 
matter and radiation. The MMR result will be valid only when the radiation temperature is much greater than the matter 
temperature, so that the recoil effect dominates the counterbalancing diffusive effect of the thermal electrons. 

In Appendix B of their paper, MMR outlined a derivation of their energy exchange formula from a Fokker-Planck 
equation that contained the Basko recoil term. In that derivation they eventually dropped the diffusion term and kept 
only the recoil terms. They clearly recognized that because of this their formula would only apply for photons of much 
higher energies than the thermal energies of the atoms. 

MMR were reluctant to include the diffusion term, because they argued that only terms of order (v /c) had been included 
and that one needed to include terms up to (v/c) 2 in order to get the reverse flow of energy from the atoms to the radiation. 
However, the present approach using detailed balance automatically provides the correct description of the reverse flow, 
as evidenced by the correct thermal equilibrium. 

6.2. Comparison with CM 

Another approach to energy exchange was given by Chen & Miralda-Escude (2004, hereafter CM), who used the F-P 
equation with Basko recoil term (14) to find the energy exchange rate in a particular cosmological context. By using this 
equation, CM assumed (as MMR did not) that this equation could be relied upon to describe the reverse flow of energy. 

CM used numerical solutions to equation (14) to evaluate the energy transfer and did not derive a general formula such 
as our equation (35), but we may do so now. In a sense the present derivation is merely extending the derivation outlined 
by MMR in their Appendix B, except that here the diffusion term is included. 

Multiplication of equation (14) by hv and integrating over v gives, after integration by parts and other straightforward 
manipulations, 

~dJ hJ] 



dU AirhkT f , , . 

■ / ^o a i» 



dt m a 



dv + kT 



dv. (37) 



In order to compare this with our result (32), we define a new "radiation temperature" Tqm by the relation 

2cm = -~. (38) 

This differs from the definition (30) by ignoring stimulated emission and also the factor v 2 in the relation between n and 
J. It is equivalent to assuming that a thermal radiation field is proportional simply to exp(~hv/kT). 

Under the same assumptions as made for equation (35), we obtain the following energy exchange formula, 

f-'G&M 1 -^- (CM) - 

We shall call this the "CM formula," because it is a natural consequence of their approach, even though they did not 
derive it explicitly. The CM formula differs from for our equation (35) by the replacement of Pro by Pcmo (and the fact 
that P does not include stimulated scatterings). 

The most noteworthy feature of the CM formula is that it shares with (35) the desirable property of taking into 
account the relationship between matter and radiation temperature (although with an slightly approximate radiation 
temperature). In particular, the approach of CM does take into account of the suppression of energy transfer for lines 
of high optical depth, in accordance with the Wouthuysen-Field effect. This clearly explains why CM found much lower 
energy exchange rates than MMR. 

Let us make some estimate of the difference between the radiation temperatures Pr and Tqm ■ Stimulated emission is 
usually fairly negligible for resonance line transfer, so let us ignore it for the moment. Then one can easily show that, 

TCM = 1-2^/hv (40) 

As an practical example, in problems involving hydrogen Lyman alpha, gas temperatures are typically below ~ 10 4 
K, since anything much higher leads to rapid ionization of hydrogen. The equivalent temperature for the Lyman alpha 
transition is ~ 10 5 K, so by equation (40) the two radiation temperatures Pr and Pcm may differ by up to ~ 20%. 
However, for lower radiation temperatures, they will differ proportionally less, and the CM radiation temperature Pcm 
may be a reasonable approximation to the more accurate Pr. For very sensitive applications, however, it would be prudent 
to use the accurate form. 



7. DISCUSSION 

A new form for the Fokker-Planck equation desribing resonant line scattering has been derived that incorporates detailed 
balance and stimulated scattering. When put into the form (27) the new equation shows a surprising unification with the 
Kompaneets equation. 

The new Fokker-Planck equation has been used to derive expressions for the rate of energy transfer between the radiation 
and the thermal gas, incorporating both recoil and the proper reverse flow of energy from the gas to the radiation. We 
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have compared our energy exchange result to those used by Madau, Meikson, & Rees (1997) and Chen & Miralda-Escude 
(2004). 

In their paper MMR noted that their energy exchange formula would be limited to cases where the radiation temperature 
greatly exceeds the gas temperature. We have verified this explicly by deriving the general dependence on the radiation and 
gas temperatures, as seen, e.g., in the factor (1 — T/Tro in our equation (35). This factor is thcrmodynamically required, 
and also will be absolutely crucial in cases of resonance line scattering of high optical depth, where the Wouthuyscn-Ficld 
effect can play a pivotal role. 

We have also derived an approximate explicit energy transfer formula closely following the work of CM, which is 
based on the Fokker-Planck equation (14). The resulting energy exchange formula (39) (which we call the CM formula) 
does approximately incorporate detailed balance and does approximately account for the relative radiation and gas 
temperatures. The CM formula should give acceptable results in most cases of interest. This demonstrates that the CM 
approach to energy transfer is clearly superior to the more limit formula of MMR. 

However, although the Fokker-Planck equation (14) and the associated CM energy exchange formula would seem to 
provide adequate results for most situations, there may be other reasons for preferring our new formulation of the Fokker- 
Planck equation (27). From a purely theoretical point of view it treats the thermodynamics in a more precise manner, and 
has the esthetically pleasing property of being unified in form with the Kompaneets equation. There may also be practical 
reasons for the using our more accurate formulation. For example, in determining the spin temperature of the 21cm line 
in cosmological contexts, one needs to consider the delicate balance of different processes, and here a thcrmodynamically 
correct description of Lyman alpha scattering could prove to be a crucial element. We plan to investigate this in the near 
future. 

The results of the present paper suggest some avenues for possible future work. The unification of the resonance line and 
electron scattering Fokker-Planck equations discussed in §5 strongly suggests that one might develop a single derivation of 
both equations, based on expansions of the transfer equations with redistribution functions that take all physical effects 
explicitly into account. Done carefully, this could resolve the somewhat unfortunate ambiguity of what the proper value 
of D{v) should be. 

One limitation of the present paper is its assumption of a homogeneous and isotropic radiation field J{v), a limitation 
shared by the Kompaneets equation. It would be desirable to derive a Fokker-Planck equation for resonance scattering 
with a spatially dependent and anisotropic radiation field, r, n), dependent on position r and direction n. This would 
be useful for cases where the incident irradiation is anisotropic, or in cases where the medium itself is inhomogencous. 

There is one way to obtain an approximate equation for anisotropic radiation fields, one that is commonly used with 
the Kompaneets equation. Neglecting stimulated emission, we write, 

ldl , Ai/£ d [ \d,J ( h 2\ 11 

This is equivalent to assuming that the extinction part of the scattering is treated with the directional intensity, but the 
emission part is treated as if the radiation field were isotropic, as in equation (20). Since scattered radiation is not too 
far from being isotropic, and becomes more so after many scatterings, this is often an acceptable approximation. We note 
that when the radiation field is homogeneous and isotropic, (41) reduces to equation (20) when integrated over all solid 
angles. 

However, equation (41) does have some definite defects. If one wanted to include stimuated scattering, the extinction 
term would have to depend on angle in some as-yet unknown way; this defect is not too serious, since stimulated scattering 
is not likely to be important for any resonance lines. A more serious defect is that the equation does not include the 
correlation of the frequency diffusion and shift on the scattering angle, an effect we know is very strong from the form of 
the angle-dependent redistribution function Rii(v, n; u', n') (Hummer 1962). 

Therefore, for theoretical reasons, as well as practical, it would be desirable to derive a Fokker-Planck equation for 
resonance scattering with a fully anisotropic radiation field. This could be done using the same correction procedure of 
the present paper, but now applied to the angle-dependent results of RD, their equation (A16). Such a derivation is 
currently under investigation. 



(41) 
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APPENDIX 



A. SYMMETRY OF THE FOKKER-PLANCK REDISTRIBUTION FUNCTION 

The redistribution function from which the RD Fokker-Planck equation (8) was derived satisfies the symmetry i?(f , v') = 
R(v' : v). However, because approximations have been made, it is not obvious that the Fokker-Planck redistribution 
function R defined implicitly through equation (8) must necessarily satisfy this same symmetry. To investigate this it is 
useful to consider the following double integral, 

G{v)R{v, v')F{v') dv dv' . (Al) 

This quantity is bilinear in the two functions F{y) and G(v), which are completely arbitrary, except that they and their 
derivatives are assumed to vanish sufficiently rapidly at the limits of integration (0, oo). 




Resonance Line Scattering 
We write the basic RD F-P result (6) for an arbitrary function F(v), 



I 



R(v,v')F(v')dv' = <p{v)F(y) + -f-jL 



D(v) 



dFjv) 
dv 



Using this for the integral over v' in (Al), we have 



G{v)R{v,v')F{v')dvdv' = I G(v) <j ip(v)F(v) + ^ 8 



2 dv 



D{v) 



dF{v) 



dv 



dv. 



Integrating by parts twice, assuming zero contributions at the limits, gives 



// 



G(v')R(v', v)F(v) dv dv' = / \<p(v)G(u) + 



Ai% d 
2 d^ 



D(v) 



dGjv) 
dv 



F(v) dv, 



(A2) 



(A3) 



(A4) 



where we have interchanged v <-> v' in the double integral on the left hand side. Interchanging F <-» G yields 

dF(v) 



dv 



D(v)- 



dv 



dv. 



(A5) 



(A6) 



J J G(v)R(v', v)F(v') dvdv' = J G{v) ^(v)F(v) + ^ 

Subtracting this equation from (A3) gives 

J J G(v) [R(v, v') - R(v', v)] F(v') dv dv' = 0. 

This relation is true for arbitrary G(v) and F(v), which implies, 

R(v,v') = R(v',v), (A7) 
that is, the RD F-P approximation maintains the exact symmetry of the redistribution function on which it was based. 

B. DERIVATION OF THE "CORRECTED" FOKKER-PLANCK EQUATION 

Since R(v, v') is symmetric, we may write for for an arbitrary function G(v), 



G{v')R{v',v)dv' = <p(v)G{y) + 



Ai/2 d 



D(v) 



dGjv) 
dv 



2 dv 

The symmetry also allows us to use equation (18) to give "corrected" versions of equations (A2) and (Bl), namely, 

d 



R c (v, v')F(v') dv' = cpF+ ^U^- 

2 dv 



dv y ' 



I 



A;y 2 8 

G(v')R c (v', v) dv' = ipG + ^U- 1 ^ 



(Bl) 

(B2) 
(B3) 



for arbitrary functions F(v) and G(v). 

Based on the preceding, we now derive a corrected F-P transfer equation. Stimulated emission can be easily incorporated 
into this derivation by starting with the full equation, 



— 157= / {[l + n(v)]R c (v,v')J(v')-[l + n(v')}R c (v',v)J(v)}dv'. 
c\ dt J 



The integrals over v' can be performed using (B2) with F = J and (B3) with G = 1 + n, which gives, 



J_9J 

cx dt 



dv 



dv 



Y 



(B4) 



(B5) 



The right hand side can be written as a frequency derivative, 

Substituting n = c 2 J/2v 2 and using equation (19), we find, after some straightforward manipulations, the "corrected" 
equation (20) in terms of the mean intensity. 



(B6) 



REFERENCES 



Adams, T.F. 1971, ApJ, 168, 575 

Basko, M.M. 1981, Astrophysics, 17, 69 

Chen, X., & Miralda-Escudc, J. 2004, ApJ, 602, 1 (CM) 

Deguchi, S., & Watson, W.D. 1985, ApJ, 290, 578 (DW) 

Field, G.B. 1958, Proc. I.R.E., 46, 240 

Field, G.B. 1959, ApJ, 129, 551 

Harrington, J.P. 1973, MNRAS, 162, 43 

Henyey, L.G. 1941, Proc. Nat. Acad. Sci., 26, 50 



Hirata, CM. 2005, preprint (astro-ph/0507102) 

Hogan, C.J., & Rees, M.J. 1979, MNRAS, 188, 791 

Hummer, D.G. 1962, MNRAS, 125,21 

Kompanccts, A.S. 1957, Sov. Phys. JETP, 4,730 

Madau, P., Meikson, A., & Rees, M. 1997, ApJ, 475, 429 (MMR) 

Peebles, P.J.E. 1968, ApJ, 153,1 

Pritchard, J.R., & Furlanctto, S.R. 2005, preprint (astro- 
ph/0508381) 



10 



Rybicki 



Rybicki, G.B., & DelFAnonio, LP. 1994, ApJ, 427, 603 (RD) 
Rybicki, G.B., & Lightman, A. P. 1979, Radiative Processes in 

Astrophysics (New York: Wiley) 
Tozzi, P., Madau, P., Meiksin, A., & Rees, M.J. 2000, ApJ, 528, 597 
Unno, W. 1952, PASJ, 4, 100 



Unno, W. 1955, PASJ, 7, 81 
Wouthuysen, S.A. 1952, AJ, 57, 31 

Zcl'dovich, Ya.,B., Kurt, V.G., & Sunyacv, R.A. 1969, Soviet Phys.- 
JETP, 28, 146 



